Phase Sensitivity of a Mach-Zehnder Interferometer 
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The best performance of a Mach-Zehnder interferometer is achieved with the input state ]A^t/2-|- 
1)|7Vt/2 - 1) + \Nt/2 - 1)\Nt/2 + 1), being Nt the total number of atoms/photons. This gives: 
i) a phase-shift error confidence 6*68% = 2.67/Nt with ii) a single interferometric measurement. 
Different input quantum states can achieve the Heisenberg scafing ~ 1/Nt but with higher prefactors 
and at the price of a statistical analysis of two or more independent measurements. 



Introduction. The central goal of interferometry is to 
estimate phase shifts with the highest sensitivity given 
a finite resource of atoms/photons. There are several 
possible interferometric configurations, which are gener- 
ally reducible to the "drosophila" Mach-Zehnder (MZ), 
see Fig. ([Pi . The MZ transforms a two-mode input state 

as \ip)out = e~'^''y^\ip)inp [ll, where the generator of the 
phase translation 9 is the y angular momentum com- 
ponent of the three-dimensional rotation group Jx — 
{a^b + aSt)/2, Jy = (a+S - ab'')/2i, J, = (a^a - StS)/2 
(a, b are bosonic annihilation operators for the two in- 
put quantum modes). How precisely can the unknown 
phase shift 9 be estimated from the output |7/')otii ? It 
is well known that, if one of the two input arms of the 
interferometer is fed with vacuum, the phase measure- 
ment uncertainty is bounded by the standard quantum 
limit (or "shot noise") A9 > 1/^/Nt, where Nt is the 
total/average number of atoms in the input state. The 
first breakthrough came in the early '80, when Caves 
showed that it is possible to reach sub shot-noise sen- 
sitivities after squeezing the vacuum of the unused port 
of the interferometer [3| . Sub shot-noise was experimen- 
tally demonstrated few years later Q- This inaugurated 
a large body of literature proposing new states and opti- 
mal performances 

d mm IE la] to push the sensitivity 
up to the Heisenberg limit ~ 1/Nt |i|. Paradoxically, 
however, fewer efforts have been devoted to the correct 
inference of phase measurement uncertainties. This fact 
went (and it is still going) almost unnoticed, despite of 
the seminal works of Helstrom and Holevo [T^ and Lane, 
Braunstein and Caves Apart from a few noticeable 
exceptions, in the current literature the phase measure- 
ment uncertainty is still retrieved from a simplified er- 
ror propagation theory, A6 — (A^)/|^^|, where (A) 
is the average of a generic phase dependent observable 
A{9) and A A = (A^) — (A)^ is the mean square fluctu- 
ation. Unfortunately, this simple approach can give cor- 
rect estimates only when the probability distributions are 
Gaussian |l2l | , which is seldom the case in quantum inter- 
ferometry. Moreover, optimized phase estimations often 
require the distribution of the available finite resource of 
Nt particles among p independent interferometric experi- 
ments, each with N = Nt/p particles. This optimization 
cannot be analyzed with the linear error propagation for- 
mula, which only considers the dependence on the num- 



ber of particles in each experiment A'' and not on the 
total number of particles Nt- As a matter of fact, even 
if the Mach-Zehnder has been extensively studied, there 
is not published literature about the exact dependence 
of phase uncertainties on the total (or mean) number of 
particles of any input state. 

Matter waves Mach-Zehnder interferometers have been 
realized with cold atoms 13] and electrons J.4j. Evi- 
dences of non-classical light properties after a beam split- 
ter have been experimentally showed in |l5j| . Entan gled 
"NOON" states have been created with few photons [13 
and ions |l7j with applications in quantum lithography 
[l8j |. Current efforts are devoted to the creation of the MZ 
with dilute Bose-Einstein condensates (BEC). EEC offers 
almost unique possibilities to create the large intensity, 
highly squeezed waves needed for Heisenberg interferom- 
etry. Recently, double-slit as well as Michelson inter- 
ferometers z(f\ have been experimentally demonstrated. 
By trapping BEC in two wells and in periodic potentials, 
the beam splitters and mirrors are replaced by the dy- 
namical tailoring of interwell barriers and magnetic wells. 
Fock and number squeezed condensates have been real- 
ized in 21]. 

In the following, we study a MZ interferometer fed by 

\tp)i7ip = -^{Ij + m)a\j -m)b + \j -m)a\j + m)b), (1) 

where N = 2j and 2m are the total and the relative num- 
ber of particles at the a and b input ports, respectively, 
cfr. Fig. (Q). We rigorously calculate error confidences 
within a Bayes framework We first consider the 

twin-Fock {m = 0) introduced by Holland and Burnett 
0. We then study a "NOON" state (m = j), showing 
that it cannot overcome shot-noise sensitivity. We will 
finally demonstrate that the best Mach-Zehnder perfor- 
mance is obtained when m = 1. This state can be created 
by constructing a twin-Fock state followed by the mea- 
surement of one particle p^ ]. 

Bayes analysis. Quantum Mechanics provides the con- 
ditional probability P{^\j,9) to measure, at the output 
ports, the relative number of atoms fj, = {Ni — N2)/2, 
given 2j — Ni + N2 particles and the unknown phase 
shift 9. With the general state Eq.J^I we have 



P(/x|j,0) = |(j,^|V)o„t| 



^m'. (2) 
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FIG. 1: Schematic representation of the Mach-Zehnder inter- 
ferometer. Atoms/photons enter the a and 6 input ports, mix 
and recombine in the beam splitters and are finally detected 
in Dl and D2. The phase shift is inferred from the number 
of atoms/photons measured in each output port. 



where |j, /J.) = \j + ^j^)Dl\j — fJ-)D2, being Dl and D2 the 
detectors at the output ports (see Fig. Q), and d^ jj{9) 
being the angular momentum matrices |23| |. The goal is 
to estimate 9 after detecting a certain value fi. Accord- 
ing to the Bayes theorem, P(9\j, fi)P{iJ,) = P{n\j, 9)P{9), 
where P{9) and take into account any further a pri- 
ori information about the real value of the phase shift and 
the output measurement ^. We assume complete igno- 
rance, so P{9) = H(7r/2 - 16*1), being H(a;) the Heaviside 
step function, while P{^) is fixed by the normalization. 
We notice that, with the state Eq.lQJ as input of the 
MZ, it is possible to estimate only the absolute value of 
phase shifts in the interval [0,7r/2]. In most cases, the 
best interferometric performance requires the statistical 
analysis of several independent measurements. Once ob- 
tained the results /ii.../ip, the phase probability distribu- 
tion reads Pp((/)|j, Mi.../ip) oc Hfc^i P{^\ii Mfe)- Averaging 
over all the possible p-uple fii...fip we find the phase dis- 
tribution for fixed 9 and j: 

{tJ,i..4ip} 

where Pp{fii...fj.p\j,9) = PifJ-k\j,9) is provided by 

quantum mechanics. As phase estimator, we choose the 
value of the phase 4ip{j,9) corresponding to the maxi- 
mum of the probability distribution Eq.Q. The phase 
uncertainty is estimated as confidence, namely the 7- 
probability that the real value of the phase shift is within 

a given interval c-y around in 0: 7 = J^^^'' d0-Pp((/)|j, 9). 

In the following we consider the case of a null phase 
shift 9 = 0. The probability distribution Eq.Q reduces 
to P{fi\j,0) = |<5^,m -I- (5^ _™p/2, where 6a,i3 is the Kro- 
necker delta. Therefore, we have two possible results 
for the measurement of the relative number of particles: 
/i = +m, and /i = —m. Because of the symmetry proper- 
ties of the MZ interferometer and of the state Eq. , we 
have that P{<j)\j,m) = P{(j)\j,—m), and the sum Eq.© 

reduces to Pp(0|j,O) « |rf^™^„,(0) + When 
9 the probability of obtaining values of /i 7^ ±m in- 
creases. In this case we might need several measurements 
in order to explore all the relevant probability distribu- 




FIG. 2: A,B: Probability distribution of the m — twin-Fock 
state with A) A'' = 40 particles and ap = 1 measurement, (A); 
A'' = 20 particles and p — 2 measurements, (B). By combin- 
ing the independent measurements it is possible to strongly 
reduce the weight of the tail of the distribution with respect 
to the central peak. C,D: Amplitude dj j^(<^) (red Une), and 

di,_i(<^) (blue line), for = 40 particles {j = 20) (C). 
The tails of the amplitudes oscillate out of phase and interfere 
destructively, so to enhance the central peak of the probability 
distribution P(<?!)|i,0) ~ |d{_i((/>) + d{ _-^{(j))\\ (D). 



tion, so the sensitivity of the interferometer for a given 
number of particles is, in general, reduced p3 |. 

Mach-Zehnder interferometer with twin-Fock States. 
We first consider the input state Itp) inp = \j)a\j)b {m = 
in Eq.^). 

The probability of measuring a relative particle num- 
ber fj, at the output port, given a phase shift 9, is |25j : 

P{fi\j,9) = g^[P/(cos6')]^ where P/(x) are the As- 
sociated Legendre Polynomials. When 9 = the proba- 
bility P(/i|j, 9) reduces to a single peak centered in n = 0. 
Therefore, repeating the experiment p times we have 
Ppi't'ljiO) ~ [P°(cos0)]^P. With a single measurement 

1 /3 

{p = 1), the 68%-confidence scales as ^ l/A^y while, 
as discussed in |2^, the mean square phase fluctuation 
scales as cr ~ 1/lniVj-. In both cases the phase uncer- 
tainty is worse than shot noise This problem is 
generally overlooked when the sensitivity is calculated 
from the simplified error propagation formula, or looking 
at the width of the central peak of the Legendre poly- 
nomials, which indeed shrink as 1/Nt- The large tails 
of the distribution contain most of the probability and 
cannot be ignored, see Fig.|2lA). This clarifies the need 
of considering several independent measurements, which 
increases the relative weight of the central peak with re- 
spect to the tails, see Fig. (|2lB). As discussed in |26j . 
the mean square phase fiuctuation is minimized combin- 
ing p = A measurements. 

In Fig. |(3lA) we study different values of confidence 
for the twin-Fock state distribution as a function of the 
number p of independent measurements and a fixed to- 
tal number of atoms Nx- We see that the bigger the 
confidence, the larger is the number of experiments, p, 
necessary to reach the minimum. It is important to no- 
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FIG. 3: Confidence interval as a function of the number of 
independent measurements p for A) the twin-Fock m = and 
B) the twin m — 1 states. Here A'r = 2000 particles. In both 
plots, from the bottom line to the top one, the black line is 
the confidence at 38.29%, the blue at 68.27%; the green at 
95.45%, the red at 99.73%, the yellow at 99.994%, and the 
brown at 99.99994%. These correspond to cr/2, cr, 2cr, 3cr, 4cr 
and 5(T when the distribution is Gaussian. 

tice that the number of measurements needed to reach 
the minimum for a given confidence does not depend on 
the total number of atoms Nt- The main results of our 
analysis can be summarized as follows: i) with a single 
measurement the MZ interferometer gives a sensitivity 
worse than the shot noise, ii) the highest sensitivity is 
reached with two independent measurements, for the 68% 
confidence, and iii) with three measurements for the 95% 
confidence: 

2.915 „ 6.654 

(4) 

The number of measurements needed to reach higher con- 
fidences is shown in Fig.||2lA). In Fig. (01 A) we summa- 
rize the scaling of the phase uncertainties reached with 
the twin-Fock state, as a function of Nt- 

How does the sensitivity scale with the number of mea- 
surements, p, given a fixed number of atoms, N , for each 
experiment? The Cramer-Rao theorem does not allow a 
better scaling than \/^fp, and we have verified numeri- 
cally that the limit is reached when p > 4, obtaining the 
mean square fluctuation « The 1/p scaling 

claimed in was probably an artifact of low statistics. 

Mach-Zehnder interferometer with NOON state. We 
now consider the state with m — j. If 6* = 
the only possible outcomes are — ±j. According 
to the Bayes theorem, we obtain the probability dis- 
tribution Pp(0|j, 0) w [cos2J(0/2) sin2J((^/2)]2p. For 
large j we have that cos^^(0/2) ~ exp[— (j/4)(^^], while 
sin^''(0/2) ~ exp[-(j/4)(0 - n)^]. Therefore, the phase 
probability distribution is a Gaussian with width scal- 
ing as cr = V2/y/T7^ (see Fig.lglB)). The NOON state 
with a Mach-Zehnder interferometer, assuming a com- 
plete "a priori" ignorance of the real phase shift, achieves 
a shot noise sensitivity. A previous analysis ^|, di- 



FIG. 4: Phase uncertainty as a function of Nt- A) Twin- 
Fock state: (from the top fine to the bottom one) black line, a 
for p = 1 ((T~ 1/ln Nt ) ; brown line, 68% confidence for p = 1 
{Ae = 0.8/7Vy'^^); yellow line, a for p = 2 (a « 1/N^^'^); violet 
line, a for p — 4 {a = 4/Nt); green line: 68% confidence for 
p = 2 (Af = 2.91/iVT). B) Phase sensitivity of the NOON 
state (blue line), a — ^2/Nt, independent of p, and of the 
twin m = 1 state (red line), AS = 2.67/Nt at p = 1. Points 
are numerical data. 

rectly addressing current experiments |l6l [T7| |. showed 
that NOON states and balanced homodyne detection can 
provide phase sensitivity scaling as l/N^^- 

Mach-Zehnder interferometer with a twin m = I state. 
We have previously shown how the large tails of the 
twin-Fock probability prevent the sensitivity to scale as 
~ 1/Nt. To kill these tails, and enhance the central 
peak of the distribution, we have to statistically ana- 
lyze several independent measurements. Here we fol- 
low a different strategy: is it possible to cancel out the 
tails by destructive interference? The answer is posi- 
tive. The tails of the probability distribution P{(f>\j, 0) ^ 
Mm,m('/') + ^m,-m('^)P obtained with the input state ^ 
cancel out when m is odd. This special feature, not 
present in the twin-Fock state, is a consequence of the 
interference between the two matrices in P((j)\j,0), and 
it can be easily recognized by expanding the Jacobi Poly- 
nomials 'Slj . Since the cancellation has a precision of the 
order of 0[{j — m)~^] we can conclude that, among this 
class of states, the case m = 1 offers the highest per- 
formance. States with larger odd values of m do not 
benefit of an equally efficient destructive interference. 
With semi-integer values of m, as for the Yurke state 
(to = 1/2) 4,], there is no cancellation at all, so we ex- 
pect that a 1 /Nt scaling can be reached only with p > 1 
as in the twin-Fock case. In Fig.||21C) we plot the am- 
plitude d{ (red line), and d{ _^{(j)) (blue line), for 

= 40 particles {j — 20). Outside the central region the 
two functions oscillate out of phase and, when summed 
to calculate P((/)|j, 0), interfere destructively. The phase 
probability distribution is shown in Fig.l2jD). The com- 
parison with the Fock state distribution (Fig. ISA)), for 
the same total number of particles, highlights the strong 
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reduction of the tails. As a consequence, for a given con- 
fidence, we reach the minimum with a smaller number 
of measurements with respect to the Fock state case, see 
Fig. |j2lB). The sensitivity of the twin-rn states is: 

2.67 , 5.376 

(5) 

The scaling of A0g8% as a function of Nt is shown in 
Fig.(@lB). The main advantage is, apart from the smaller 
prefactors, that the 1/Nt scaling is reached with a sin- 
gle — 1) measurement for both the 68% and the 95% 
confidence, as compared with the p = 2 and p = 3 mea- 
surements requested with Fock states, respectively. For 
higher confidences, the required number of measurements 
is shown in Fig. 06) [Hill. 



Conclusions. There are not known protocols for the 
direct measurement of phase shifts. Phases can only be 
inferred after the measurement of a different (phase de- 
pendent) observable. What is the highest achievable sen- 
sitivity? Within a rigorous Bayesian analysis we find that 

\i')^np - l/V2(|iVT/2 + l)|iVT/2-l) + |iVT/2-l>|A^T/2 + 

1)), gives the highest sensitivity Cg8% — 2.67/Nt, with 
a single interferometric measurement. Different input 
quantum states can achieve the Heisenberg scaling ~ 
1/Nt but with higher prefactors and at the price of a 
statistical analysis of two or more independent Mach- 
Zehnder measurements. 
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The possibility to estimate a phase with a single mea- 
surement has some important advantages. Indeed, when 
the external perturbation has a short life-time, the need 
to perform several independent measurements might re- 
quire the building of an equal number of interferometers. 
We have extended our analysis to consider classes of 
states different from Eq.Q. We have always found that 
the twin m = 1 state provides the best performance, 
and, on the ground of some general arguments, we be- 
lieve that this remains true with any possible two-mode 
input quantum state. 



